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Abstract 

A method of computing Case's eigenvalues is proposed. The eigenval- 
ues are obtained as eigenvalues of a tridiagonal matrix. 



I ■ 1 Introduction 

j^ ' In this paper we consider eigenvalues in Case's method [1] of solving the radiative 

transport equation. For isotropic scattering, discrete eigenvalues v of Case's 
method are given by solutions to 1 — ci^tanh" (l/i^) = 0, where c is a constant 
[9]. In the general case of anisotropic scattering, the largest eigenvalue can 
be found by Holte's expansion [21 [31 IH [S]. The behavior of roots has been 
intensively studied [S1[71[S]. However, we have many discrete eigenvalues when 
^^ ■ scattering is highly anisotropic such as the Henyey-Greenstein model [5]. Then, 

fN| I finding discrete eigenvalues as roots is not easy. Here we compute eigenvalues 

VO ' making use of a tridiagonal matrix, which we refer to as Markel's matrix i?™ 

tH- ' [lOl [11] . All eigenvalues are obtained by one diagonalization. 

o' 
m 

2 Case's method 

Let u{x,s) be the specific intensity at position x e R in direction s e §^. Let 
^^ , /x be the cosine of the polar angle with respect to the positive z axis, and (p be 

the azimuthal angle. The parameter c > 0. The (time-independent) radiative 
transport equation is given by 

fi—-u(x,s)+u(x,s) — c /(s • s')u(a;,s')ds', (1) 

ox Js2 

We give the phase function /(§ • s') as a polynomial of order N: 

N I 

/(s.§')=E E fiYim{s)Y;„,{s'). (2) 

1=0 m=-l 
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In the case of isotropic scattering, N = and /o = 1. In the Henyey-Greenstein 
model [8], TV = cx) and /; = g', where g G [—1, 1] is a constant. 
We seek particular solutions of ([T]) of the form 

u(a;,s) = $™(s)e--/^ (3) 

where [12] 

$™(s) = ^™(i.,/x)(l-/i2)l™l/%W^ J 0™(^.,^)(l-^2yH^^^^^ (4) 

We obtain 

(l-a$™(s) = c/ /(s-sO^rCSOds'. (5) 

We define 

ai^l~cfiQ{N~l), (6) 

where the step function &{■) is defined as 9(x) = 1 for x > and — for a; < 0. 
We introduce polynomials /i["(i/) (0 < / and |m| < TV) as [HITS] 



'H=//'"(-,A^)(l-A^^)""'^'(-l) 



We then have the recurrence relation, 

i^{2l + l)(Jih]"{i') -{l-m+ l)h]'^^(i^) -{1 + m)hY'_^{v) ^ 0, (8) 

with 

>Al,,\ _ ^r)|_| iMi _ (2|w|). |„| / x _ /„, I , -,N_, ,|m|. 



/.j:|H = (2H-l)!! = ^_L_i^, C|+i(-) = (2H + 1Wh/z|:j(.). (9) 

We also have 

/,-|™l(.) = (-l)l™l|_H|;,l™l(,). (10) 

Let us define 

^ (]' — V 

7"(z^,m) = (-1)" E /r(2r + l)^— ^P,?'(A.)(l-M')"""'^'/irM- (11) 
i'=\m\ ^ -\-m). 

Singular eigenfunctions (j)"^ {1^,11) are obtained as 

^-(:., /i) = ^T'^;^^^^^ + A'"H(1 - ^2)-l™l^(^ „ ^), (12) 

2 V — fJL 

where 






The continuous spectrum lies on [—1,1], and discrete eigenvalues v'^ {j = 
0, . . . , M — 1) are computed as roots of A'"(z), z S C: 

A'"(z) = 0, (14) 

where 

A'"(z) = i-£f /■' 2:!^(fii^(i_ ^2)1-1 ^^_ (15) 

2 J _i z — fi 

The number M of discrete eigenvalues depends on \m\ and we have [HI [T^] 

M <N - \m\ + 1. 

3 Expansion by spherical harmonics 

Finding roots by (|14l) is not easy when scattering is highly anisotropic {N is 
large). Here we compute eigenvalues making use of Markel's matrix B"^ [TOl[TT| . 
which was first introduced in the method of rotated reference frames. Let us 
expand singular eigenfunctions as [TH [TB] 

oo 

$™(s) = Y. <^T{^)yim{^- (16) 

/ — I m 

The coefficients c"^{y) are calculated using (IT2|) . We multiply y;^(s) on both 
sides of ([5]) and integrate them with respect to s. 

/ (l-a$™(s)y,:;,(s)ds = c/ / /(s.s')$™(s')l^;,(s)dsds'. (17) 

We obtain 

cr(^) - - E f / A^>1'™(§)i^™(S)ds) cl?(;.) = c/,e(iV - Ocr(^)- (18) 

Hence we arrive at an eigenproblem: 

B''^\^"^{v)) =v\%Ij'^{v)), (19) 

where 



sj? = ^^ / A*>^p™(s)r,:;,(s)ds 



'^'"^^ -^^■^-^^ IJXII^':^ ^^'..1. (20) 



(4^2 - l)aza,_i -'-^ V (4(^ + 1)' - l)'^i+i'Ti 



(^IV'"(^)> = ^===^/^c^(^), (21) 

V^ (^) 
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Figure 1: Isotropic scattering with iV = 0. We set c = 0.9. Eigenvalues v for 
different Nb are shown. 



where Z™{iy) is the normahzation constant to guarantee {ip"^{i')\ip™{i^)) — 1. 
Therefore, we can compute v as eigenvalues of B"^. If < c < 1, z^ are real 
values because B™ is a real symmetric matrix. If i' is an eigenvalue, then —ly is 
also an eigenvalue [10]. Note that v depends on m, and |i^™| > 1. In particular, 
discrete eigenvalues i^J^ are obtained using B™. 

To numerically calculate discrete eigenvalues, we can diagonalize an Nb x Nb 
matrix B™ with Nb sufficiently large that discrete eigenvalues are obtained 
with desired accuracy. Figures [T] and [2] show eigenvalues for different Nb in the 
cases of isotropic scattering and the Henyey-Greenstein model with g — 0.9 and 
m = 0, respectively. In Fig. |3l eigenvalues for g = 0.9 are plotted for different 
m. As is shown in [IB] , the value of the largest eigenvalue sharply drops as m 
increases. 



4 Conclusions 



In this paper, we have proposed a novel method of computing Case's discrete 
eigenvalues using a tridiagonal matrix i3™. Discrete eigenvalues are computed 
by one diagonalization. In Figs. [T] and [2l we see that discrete eigenvalues are 
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Figure 2: The Henyey-Greenstein model with N ^ 9, g ~ 0.9. We set c = 0.9. 
Eigenvalues ly for different A'^s are shown for m = 0. 
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Figure 3: The Henyey-Greenstein model with N ~ 9, g — 0.9. We set c = 0.9. 
Eigenvalues ly for different m are shown for iV^ = 501. 



quite stable as Nb increases; relatively small matrices give accurate discrete 
eigenvalues. For isotropic scattering A^ = 0, there are two discrete eigenvalues 
ztj/Q. If c = 0.9, then we find vq — 1.9032 from the transcendental equation 
1 — ci'o tanh~^(l/i^o) = 0. On the other hand, the largest eigenvalue of B^ is 
obtained as 1.8257, 1.9027, 1.9032, 1.9032 for Nb = 1, 3, 5, and 501, respectively. 
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